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We present a measurement of the top-quark mass, M top , in the all-hadronic decay channel ti — > 
W + bW~b — > qiq2b<?3</4&- The analysis is performed using 310 pb -1 of ^/s=1.96 TeV pp collisions 
collected with the CDF II detector using a multi-jet trigger. The mass measurement is based on an 
event-by-event likelihood which depends on both the sample purity and the value of the top-quark 
mass, using 90 possible jet-to-parton assignments in the six-jet final state. The joint likelihood of 
290 selected events yields a value of M top =177.1 ± 4.9 (stat.) ± 4.7 (syst.) GeV/c 2 . 

PACS numbers: 14.65.Ha, 12.15.Ff 
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The mass of the top quark, M top , is an important free 
parameter of the standard model (SM) and is of the order 
of the electroweak symmetry breaking scale. Since vir- 
tual top quarks are involved in higher-order electroweak 
processes, by measuring the top-quark mass one can con- 
strain the mass of the SM Higgs boson [l[ and particles 
predicted in extensions of the SM Q. At the Tevatron, 
ti pairs are produced by QCD processes and decay ac- 
cording to the CKM matrix [3] with a branching ratio 
of almost 100% into a W boson and a b quark. The fi- 
nal state of the event is then defined by the decay of the 
W bosons. All-hadronic ti events, where both of the W 
bosons decay into quarks, have a multi-jet final state and 
no missing energy due to neutrinos. The top-quark mass 
measurement in this decay channel is motivated by the 
large branching fraction (w44%) and the complete recon- 
struction of both top and antitop quarks, relying only on 
hadronic jets. It is the first top-quark mass measurement 
in this channel using Tevatron Run II data, and serves as 
a consistency check of the measurements in the two other 
ti decay modes 0,0]. The major experimental challenge 
is the presence of a large amount of multi-jet background 
events from QCD processes, which dominates ti produc- 
tion by three orders of magnitude, even after applying a 
trigger dedicated to events with a multi-jet final state. 

In this Letter, we present the measurement of the top- 
quark mass in the all-hadronic decay channel using a 
sample of ti decays corresponding to 310 pb _1 of proton- 
antiproton collisions at y/s—1.96 TeV, collected using the 
CDF II detector between February, 2002, and August, 
2004. We measure the mass using the ideogram method, 
which was used in the DELPHI experiment at LEP for 
the W boson mass measurement |6j . The method is based 
on an event-by-event likelihood that reconstructs the top 
and antitop quarks using all 90 assignments of jets to 
quarks from W boson decays and b quarks (jet combina- 
tions) from the ti decay. Considering all jet combinations 
enhances the amount of top-quark mass information we 
can extract from each event. The 90 jet combinations 
arise from the 6!/8 permutations of the six jets, where 
the two possible top antitop permutations and the two 
possible permutations of each of the two quark pairs com- 
ing from the W decay are treated equally. 

The sample purity, V, defined as the fraction of all- 
hadronic ti events contained in the selected data sample, 
is a free parameter of the likelihood together with the 
value of M top . The reason for measuring both M top and 
V simultaneously is that the cross sections of the QCD 
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backgrounds are not well known. Even though the back- 
ground processes do not contain any the top-quark mass 
information, the sensitivity of the likelihood on M top does 
depend on the total amount of QCD background events. 
This is reflected both in the statistical variance of M top , 
and in corrections applied to M top which take into ac- 
count the presence of background events. 

The CDF II detector Q is a general-purpose charged 
and neutral particle detector designed to study pp colli- 
sions at the Fermilab Tevatron. It consists of an eight- 
layer silicon microstrip detector array and a cylindrical 
drift chamber contained inside a 1.4 -T solenoid magnet, 
surrounded by electromagnetic and hadronic sampling 
calorimeters with a geometrical acceptance up to a pseu- 
dorapidity of = 3.6 (|. Muon chambers are located 
outside the calorimeters, allowing the reconstruction of 
track segments for penetrating particles. 

The dataset used for this measurement is selected with 
a multi-jet trigger that relies solely on calorimetry. This 
trigger requires at least four calorimeter clusters with a 
transverse energy Ei > 15 GeV and the scalar sum of the 
transverse jet energies to exceed 125 GeV. To improve the 
signal-to-background ratio, we impose offline kinematic 
requirements based on the scalar sum of the transverse 
jet energies and event shape observables. A full descrip- 
tion of the trigger and kinematic selection is given in an 
earlier publication [i| . For the final top-quark mass mea- 
surement we consider only events with exactly six jets, 
each with transverse energy Et > 15 GeV and pseudo- 
rapidity 1 77 1 < 2. Jets are identified as clusters of energy 
deposits in the calorimeter segments (towers) which fall 
within a cone radius of AR = ^Acj) 2 + A77 2 < 0.4 
The raw jet energy scale (JES) is the result of a multi-step 
correction procedure to convert the measured transverse 
jet energy to the expected transverse energy of the par- 
ton corresponding to the jet. The corrections, assessed 
using data and a simulation of the CDF detector, include 
corrections for the response inhomogeneity in 77, contri- 
butions from multiple interactions, the non-linearity of 
the calorimeter response, the underlying event, and the 
energy flow out of the jet cone [101 ] . We apply addi- 
tional jet-parton energy corrections, specific to ti events, 
which are parameterized independently for 6-quark jets 
and light-flavor jets 

In order to further reduce the QCD background, we 
apply 6-quark identification to each event passing the 
kinematic selection. At least one jet must have a recon- 
structed displaced vertex (6-tag), consistent with a long- 
lived bottom hadron [I2I ] . This displaced vertex informa- 
tion is not used to reduce the number of jet combinations 
in the event likelihood. Instead, we derive a weight factor 
using the CDF jet probability algorithm [l]| which takes 
into account the probability that a jet originates from a 
b quark, as explained below. The final data sample con- 
tains 290 events. Assuming a theoretical cross section of 
6.1 pb [14] for a top-quark mass of 178 GeV/c 2 15j, the 
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expected signal-to-background ratio is about one-to-five. 

By applying energy and momentum conservation we fit 
the momenta of the jets and the two top-quark masses, 
m, 1,2 , and determine their estimated uncertainties, cr*' 2 
using a x 2 minimization. This kinematic fit is applied to 
each of the 90 jet combinations, i, where we constrain the 
masses of the two W bosons to a world average value of 
80.4 GeV/c 2 Q within their natural widths. In or- 
der to achieve a better background reduction, the two re- 
constructed top-quark masses, m|' 2 , are not constrained 
to be equal. The m.' distributions thus obtained peak 
around 160 GeV/c 2 and range from 100-300 GeV/c 2 for 
the whole data sample, while the uncertainties cr*' 2 are 
typically of the order of 15 GeV/c 2 and range from 5- 
30 GeV/c 2 . We measure M top using an event likelihood 
where we weigh each combination i with a factor «;,, ex- 
pressing the compatibility with both the tt decay kine- 
matics and the presence of two 6-quark jets together with 
four light-flavor jets. The weight factor wt is calculated 
as 

t«i=exp(™Xi)ti>i, (1) 

where the exponential term, calculated using the \ 2 value 
of the kinematic fit, is a measure of the compatibility of 
jet combination i with the kinematics of a decaying tt 
pair. The second factor, w\, gives the probability of two 
of the jets to originate from a b (or b) quark and four 
others to originate from light quarks (including c quarks), 
and is given by 

where we assume quarks with indices 3-6 to be light 
quarks originating from the W bosons. The probability 
pf° T ^ for a jet j to originate from a b quark (or light- 
flavor quark, q) is obtained using the CDF jet probability 
algorithm [13j. This calculates the probability Pj for a 
given jet j to originate from the primary vertex, based 
on the impact parameter information from the tracks be- 
longing to that jet. The Pj distribution for &-quark jets 
(or light-flavor jets) is described by a probability density 
function, f(b(or q)\Pj) in order to obtain 

b(or g) = /(&(or q)\Pj) ,„» 

Pj fibip^+fmr [ 1 

The event likelihood consists of two terms. The sig- 
nal likelihood term corresponds to the convolution of 
two Breit-Wigner distributions, FbwK? 7 ^ |M top ), with 
two Gaussians, G(m'j\m\,<jl), describing the experimen- 
tal resolutions, , for each of the two reconstructed 
top-quark masses, m*' 2 : 

CWo P )= n /gKKX) 



xF BW (m' 3 \M top )dm' 3 . (4) 

The Mtop-independent background likelihood term, C^ s , 
corresponds to the two-dimensional posteriori probability 
density function (histogram) of m 1 ^ 2 obtained from ALP- 
GEN [lj] Monte Carlo (MC) multi-jet QCD background. 
All MC events are passed through the CDF detector sim- 
ulation and are subjected to the same event selection cri- 
teria as the data. 

The likelihood for a given event n is derived by sum- 
ming the signal and background event likelihoods for each 
jet combination i and is calculated as a function of the 
top-quark mass, M top , and the sample purity, V: 

90 

C"(M top ,V) =J2wi[PjCf S (Mto P ) + (1 - nC<*]- (5) 

2=1 

We obtain a one-dimensional likelihood curve as a func- 
tion of Mtop by maximizing the two-dimensional joint 
likelihood, 

290 

L(M top ,V) = [] C n (M top ,V), (6) 
n=l 

with respect to the sample purity, V, for each value of 
M top . The maximum of the total likelihood of the 290 se- 
lected data events corresponds to a sample purity value 
of V — 0.21 ± 0.07(stat.), which is compatible with a 
signal-to-background ratio of about one-to-five expected 
from SM tt production. The value of M top extracted 
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FIG. 1: The correlation between the measured value ol Mtop 
and the generated top-quark mass obtained with pseudo- 
experiments using tt signal events generated with herwig and 
an admixture ol QCD background generated with ALPGEN. 
The solid line represents the linear parameterization used to 
correct Af top . The shaded band shows the effect ol a one- 
sigma variation ol the sample purity, V = 0.21 ± 0.07(stat.). 
The dashed line corresponds to 100% correlation and zero 
offset. 

from the joint likelihood fit in Eq. ^ is biased due to the 
presence of wrong jet combinations, background events, 
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FIG. 2: The one-dimensional negative log-likelihood of M top 
obtained including corrections using the parameterization 
shown in Fig. [T] 



and assignment of jets that arise from initial- (ISR) and 
final-state gluon radiation (FSR). Therefore, we correct 
the measurement of M top using a linear parameteriza- 
tion of the correlation between the generated top-quark 
mass value and the mass estimator M top in the range be- 
tween 160 GeV/c 2 and 190 GeV/c 2 , as shown in Fig. Q] 
Each mass point is obtained from pseudo-experiments 
with MC samples containing a mixture of 290 Poisson- 
fluctuated signal plus background events. We use her- 
WIG [l3| to generate the ti signal events and ALPGEN 
to model the QCD background. The tag rate matrix 
used in models the false identification of 6-quarks. 
The correction depends on the sample purity used to 
construct the pseudo-experiments as illustrated by the 
shaded band which corresponds to the one-sigma statis- 
tical uncertainty on the sample purity measured from 
data. Accordingly, we assign a systematic uncertainty 
due to this purity uncertainty. Using the same pseudo- 
experiments we verify that the measured sample purity is 
independent of the generator top-quark mass value and 
that the statistical uncertainty on M top covers 68% of the 
sample measurements. We correct the statistical uncer- 
tainty from the likelihood curve using the width of the 
pull distribution (M top — ro g en)/cM top , where m gen is the 
generated top-quark mass, M top the measured top-quark 
mass, and OM to its statistical uncertainty obtained in 
the pseudo-experiments. The correction factor, 1.17, is 
independent of the value of M top . 

The maximum of the one-dimensional top mass likeli- 
hood, after correction with the parameterization shown 
in Fig. [T] corresponds to the value M top = 177.1 ± 
4.9(stat.) GeV/c 2 . Figure. [2] shows the corresponding 
likelihood. 

The uncertainty on the JES which varies between ap- 
proximately 3% and 8% of the measured jet energy, de- 
pending on the r\ and pt of the jet [l(| , dominates the sys- 



tematic uncertainties on our measurement of M 



top ■ 



We 



determine its effect using pseudo-experiments in which 



the JES is increased (decreased) with one standard de- 
viation. Half of the difference in M top between the 
positive and negative variation of the JES amounts to 
AM top = 4.3 GeV/c 2 . 

The corrections to M top , shown in Fig. [TJ depend both 
on the background fraction and shape. The uncertainty 
due to the background fraction, varied within the un- 
certainties of the measured sample purity, is AAf top = 
1.1 GeV/c 2 for the measured value of M top . The relative 
fractions of the light-flavor QCD background and QCD 
background containing b quarks is not precisely known. 
We estimate their ratio (light-flavor/6) to be 2.8 ± 1.0 
by performing a binned maximum likelihood fit of the 
MC expectation to the data based on the jet probability 
weight distribution in Eq. ([2]). This result is consistent 
with a value of 2.7 obtained using the respective QCD 
cross sections given by the ALPGEN event generator, tak- 
ing the individual selection efficiencies into account. A 
systematic uncertainty due to the background composi- 
tion results in AM top = 0.8 GeV/c 2 by varying the back- 
ground flavor composition between only light-flavor QCD 
background events and only QCD background containing 
b quarks. 

In accordance with || , we vary two pythia [3] parton 
shower parameters, Aqcd and the ISR/FSR transverse 
momentum scale X-factor, in order to model the system- 
atic uncertainties due to the amount of ISR and FSR. The 
largest difference observed in M top between any variation 
amounts to 0.9 GeV/c 2 (0.8 GeV/c 2 ) for ISR (FSR). 

To examine the systematic effect due to uncertainties 
in the parton distribution function (PDF) of the proton 
and antiproton, we follow the approach used in [5(, re- 
sulting in AMt op = 0.8 GeV/c 2 . 

By applying the parameterized corrections on M top ob- 
tained with herwig to MC events generated with pythia 
we estimate MC modeling uncertainties. The resulting 
difference in M top equals 0.5 GeV/c 2 . 

Uncertainties in &-quark fragmentation, semi-leptonic 
branching ratios, and color-flow affect the 6-quark 
JES This results in an additional contribution to 
AM top of 0.5 GeV/c 2 . All systematic uncertainties are 
summarized in Table [T] They result in a total systematic 
uncertainty of AM top = 4.7 GeV/c 2 . 

In summary, we present a new measurement of the top- 
quark mass in the all-hadronic channel. The measured 
value is M top = 177.1 ± 4.9(stat.) ± 4.7(syst.) GeV/c 2 . 
This measurement is the first determination of the top- 
quark mass in the all-hadronic channel using Run II data 
and is twice as precise as the Run I measurements in this 
channel H, 2$, ■ We expect the statistical uncertainty of 
this measurement to be reduced to 2 GeV/c 2 after ana- 
lyzing 1 fb _1 of integrated luminosity. Our result can be 
compared with the latest CDF top-quark mass measure- 
ments in the dilepton 4] and lepton+jets [H channel, 
which are based on samples with a similar integrated 
luminosity. The measured value for the dilepton chan- 
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TABLE I: Systematic uncertainties on the measured top- 
quark mass. 



Systematic 
uncertainty source 


AM top [GeV/c 2 ] 


[3] 


JES 


4.3 




Background fraction 


1.1 


[4] 


Background shape 


0.8 




ISR 


0.9 




FSR 


0.8 


[5] 


PDF 


0.8 


MC modeling (generator) 


0.5 




6-quark JES 


0.5 


[6] 



Total 



4.7 



nel is Af top = 167.9 ± 5.2(stat.) ± 3.7(syst.) GeV/c 2 , 
and the value for the lepton+jets channel is M top = 
173.5l|;g(stat. + JES) ± 1.3(syst.) GeV/c 2 , where the 
first uncertainty includes the uncertainty on the jet en- 
ergy scale, which is measured simultaneously with M top . 
At this level of precision, the values of M top are compat- 
ible between all three decay channels. 
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